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Abstract 

We prove the existence of both local and global smooth solutions to the 
Cauchy problem in the whole space and the periodic problem in the n-dimensional 
torus for the incompressible viscoelastic system of Oldroyd-B type in the case of 
near equilibrium initial data. The results hold in both two and three dimensional 
spaces. The results and methods presented in this paper are also valid for a wide 
range of elastic complex fluids, such as magnetohydrodynamics, liquid crystals 
and mixture problems. 

1 Introduction 

Many of the rheological and hydrodynamical properties of complex fluids can be at- 
tributed to the competition between the kinetic energies and the internal elastic ener- 
gies, through the special transport properties of their respective internal elastic vari- 
ables. Moreover, any distortion of microstructures, patterns or configurations in the 
dynamical flow will involve the deformation tensor F . In contrast to the classical sim- 
ple fluids, where the internal energies can be determined by solely the determinant 
of the deformation tensor F, the internal energies of the complex fluids carry all the 
information of this tensor [19], [8]. 

In this paper we consider the following system describing incompressible viscoelas- 
tic fluids. The existence results we obtain in this paper, together with the methods, 
are valid in many related systems, such as those for general polymeric materials O [19] , 
magnetohydrodynamics (MHD) [9], liquid crystals [lOlEl], and the free interface mo- 
tion in mixture problems [38] . The entire coupled hydrodynamical system we consider 
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here contains a linear momentum equation (force balance law), the incompressibility 
and a microscopic equation specifying the special transport of the elastic variable F: 



V ■ = 0, 

vt + v■Vv + Vp = ^l^v + V■[^^F^l (1.1) 
Ft + v-VF = VvF. 

Here v{t,x) represents the velocity field of materials, p(t,x) the pressure, fi (> 0) the 
viscosity, F{t,x) the deformation tensor and W{F) the elastic energy functional. The 
third equation is simply the consequence of the chain law. It can also be regarded as 
the consistence condition of the flow trajectories obtained from the velocity field v and 
those from the deformation tensor F fiTi [25| [271 [8] . Moreover, in the right hand side 
of the momentum equation, ^^^^^ is the Piola-Kirchhoff stress tensor and ^^^^F^ is 
the Cauchy-Green tensor, both in the incompressible case. The latter is the change 
variable (from Lagrangian coordinate to Eulerian coordinate) form of the former one. 
Throughout this paper we will adopt the notations of 

Ov ■ 

iVv)ij = {VvF)i, = (yv)ikFkj, (V ■ F), = 

and summation over repeated indices will always be well understood. 

The above system is equivalent to the usual Oldroyd-B model for viscoelastic fluids 
in infinite Weissenberg number cases [I9]. On the other hand, without the viscosity 
term, it represents exactly the incompressible elasticity in Eulerian coordinate. We 
want to refer to [211 |22l |23l [261 [271 HI El E] and their references for the detailed 
derivation and physical background of the above system. 

Due to the elasticity nature of our system (also being regarded as a first step 
in understanding the dynamical properties of such systems), the study of the near 
equilibrium dynamics of the system is both relevant and very important. For this 
purpose, we will impose the following initial conditions on system (11.11) : 

F{0,x) = I + Eo{x), v{0,x) =Vo{x), x e Q. (1.2) 

where Q is the physical domain under consideration. We further assume that Eq{x) 
and Vo{x) satisfy the following constraints: 



V-Vo = 0, 
det(/ + Eq 



V-E^ = 0, ' (^-^^ 

VmEoij — VjEoim = EqijV lEoim — EoimVlEoij. 

The first three are just the consequences of the incompressibility condition ^271 [25] and 
the last one can be understood as the consistency condition for changing of variables 
between the Lagrangian and Eulerian coordinates (see Lemma [2.41 and Remark 12. 5p . 
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When i7 is a bounded domain with smooth boundary, we will choose the following 
Dirichlet boundary conditions: 

v{t, x) = 0, E{t, x) = 0, {t, x) e [0, T) X dQ. (1.4) 

Global existence of classical solutions for system (11 .ip with small initial data Eq 
and Vq, for the Cauchy problem in the whole space and the periodic problem in the 
n-dimensional torus f2 = T", will be proved in this paper. Our methods in this paper 
are independent of the space dimensions. We point out that the initial-boundary value 
problem (11.11) . with (11.21) and (11.41) can also be treated at a more lengthy procedure, 
with few more technical difficulties than the ones presented in this paper. 

There have been a long history of studies in understanding different phenomena for 
non- Newtonian fluids, such as those of Erichsen-Rivlin models [HI [36] , the high-grade 
fluid models [121 EH EH] and the Ladyzhenskaya models [18] . There is an important 
difference between the system (11.11) considered here and all the models mentioned 
above, namely, the system (11.11) is an only partially dissipative system. This brings 
extra difficulties in the usual existence results for small data global solutions. 

There also exists a vast literature in the study of compressible nonlinear elasticity 
[U [32] and nonlinear wave equations [3, [H [TJ [151 [13 [33] . The powerful techniques, 
the generalized energy methods, which involve the rotation, Lorentz and scaling invari- 
ance, were originally developed by John and Klainerman for studying the solutions to 
nonlinear wave equations [H]. The method was later generalized by Klainerman and 
Sideris to the nonrelativistic wave equations and elasticity equations with a smaller 
number of generators, with the absence of the Lorentz invariance [ITJ [32] . However, in 
the case of viscoelasticity, the presence of the viscosity term Av prevents the system 
from possessing the scaling invariant properties. Moreover, the incompressibility is in 
direct violation of the Lorentz invariant properties [Ml [35] . 

In the compressible nonlinear elasticity, the special null condition on the energy 
functional W{F) (or the nonlinear term in the nonlinear wave equations) has to be 
imposed to carry out the dispersive estimates for the classical solutions. Due to the 
presence of the viscosity term Av, no attempt has been made in this paper to establish 
the dispersive estimates or to understand the nonlinear wave interaction/cancellations 
using the null conditions in these those in [321 [31] (although they are under 

investigation). In fact, we use a kind of standard energy estimate as those used for 
the Navier-Stokes equations. The methods in this paper are the higher order energy 
estimates, which take advantages of the presence of the dissipative term Av in the 
momentum equation and do not take into account of the null conditions on the elastic 
energy function W{F). However, due to the absence of the damping mechanism in 
the transport equation of F, we have to use some special treatment which involves 
the revealing of the special physical structures of the system. Notice that the usual 
energy method [TBI BOj does not yield the small data global existence, since there is no 
dissipation on the deformation tensor F. Motivated by the basic energy law (see the 
next section) and our earlier work in 2-D cases [251 [22], we analyze the induced stress 
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term. After the usual expansion around the equihbrium, we notice that V ■ F does 
provide some weak dissipation. 

The other key ingredient in this paper is the observation that V x F is a high order 
term for initial data under our physical considerations. Formally, this is merely the 
statement that the Lagrangian partial derivatives commute. Lemma 12.41 demonstrates 
the validity of this result in the evolution dynamics of the PDE system. 

The small data global existence of the classical solutions for the incompressible vis- 
coelastic system (11. ip provides us better physical understanding of this general system. 
The proof of the theorems involves all the special coupling between the transport and 
the induced stress, the incompressibility and the near equilibrium expansion. More- 
over, the bounds for the initial data (which depends on the viscosity) may also shed 
some lights on the large Weinssenberg number problem in viscoelasticity. 

As for the related work on the existence of solutions to nonlinear elastic (without 
viscosity) systems, there are works by Sideris [32] and Agemi [1] on the global existence 
of classical small solutions to 3-D compressible elasticity under the assumption that the 
nonlinear terms satisfy the null conditions. The former utilized the generalized energy 
method together with the additional weighted estimates, while the latter's proof 
relies on the direct estimations of the fundamental solutions. The global existence for 
3-D incompressible elasticity was then proved via the incompressible limit method [3lj 
and very recently by a different method [35]. It is worth noticing that they used an 
Eulerian description of the problem, which is equivalent to that in [271 [25] . Global 
existence for the corresponding 2-D problem is still open, and the related sharpest 
results can be viewed in [3l H]. For incompressible viscoelastic fluids, Lin, Liu and 
Zhang [25] proved the global existence in 2-D case, by introducing an auxiliary vector 
field as the replacement of the transport variable F. Their procedure illustrates the 
intrinsic nature of weak dissipation of the induced stress tensor. Lei and Zhou [22] 
obtained the same results via the incompressible limit where they directly worked on 
the deformation tensor F. Recently Lei, Liu and Zhou [21] proved global existence for 
2-D small strain viscoelasticity, without assumptions on the smallness of the rotational 
part of the initial deformation tensor. Finally, after the completion of this paper, we 
became aware of the manuscript [6] which studied the similar problems as in this paper. 

The paper is organized as follows. In section 2, we review some of the basic concepts 
in mechanics. Some important properties in both fluid and elastic mechanics will also 
be presented. Section 3 is devoted to proving local existence. The proof of global 
existence is completed in section 4. In section 5, the incompressible limit is studied. The 
result may be important for the study of numerical simulations and other engineering 
applications. 

2 Basic Mechanics of Viscoelasticity 

In this section, we will explore some of the intrinsic properties of the viscoelastic 
system presented at the beginning of the paper. These properties reflect the underlying 
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physical origin of the problem and in the meantime, are essential to the proof of the 
global existence result here. 

We recall the definition of the deformation tensor F. The dynamics of any mechan- 
ical problem (under a velocity field), no matter in fluids or solids, can be described 
by the flow map, a time dependent family of orientation preserving diffeomorphisms 
x{t,X), < t < T. The material point (labelling) X in the reference configuration is 
deformed to the spatial position x{t, X) at time t, which is in the observer's coordinate. 

The velocity field v{t, x) determines the flow map, hence the whole dynamics. How- 
ever, in order to describe the changing of any configuration or patterns during such 
dynamical processes, we need to define the deformation tensor F{t,X): 

F{t,X) = —{t,X). (2.1) 

Notice that this quantity is defined in the Lagrangian material coordinate. Obviously 
it satisfies the following rule [TT] : 

dF{t,X) _ dv{t,x{t,X)) 

dt ~ dx ■ ^^-^^ 

In the Eulerian coordinate, the corresponding deformation tensor F{t, x) will be defined 
as F{t,x{t^X)) = F(t,X). The equation (12. 2p will be accordingly transformed into 
the third equation in system (II. ip through the chain rule [IHl Hi], [27] . In the context 
of the system, it can also be interpreted as the consistency of the flow maps generated 
by the velocity fleld v and deformation fleld F. 

The difference between fluids and solids lies in the fact that in fluids, the internal 
energy can be determined solely by the determinant part of F (through density) and 
in elasticity, the energy depends on the whole F. 

The incompressibility can be exactly represented as 

detF = l. (2.3) 

The usual incompressible condition V ■ f = 0, the flrst equation in (ll.ip . is the direct 
consequence of this identity. 

Since we are interested in small solutions, we deflne the usual strain tensor in the 
form of 

E = F-L (2.4) 

The following lemma is well known and appeared in [37] . It illustrates the incom- 
pressible consistence of the the system (II. ip . 



Lemma 2.1. Assume that the second equality of 1^1. 3\} is satisfied and {v,F) is the 
solution of system ( fi.i)) . Then the following is always true: 

det{I + E) = l (2.5) 

for all time t > 0. 
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Proof. Using the identity ^^r^ = (detF)F the first and third equations of fll.ip 
give the result 

( det(/ + E))^ + V ■ V( det(/ + E)) 

= det(J + E){I + E)jlVkVi{I + E)kj 
= det{I + E)V -v = 0. 

Thus, the proof of Lemma 12.11 is completed. □ 

The following lemma played a crucial rule in our earlier work [251 [2? J . It provides 
the third equation in (11.11) with a div-curl structure of compensate compactness [27], 
as that in the vorticity equation of 3-D incompressible Euler equations. 



Lemma 2.2. Assume that the third equality of U.!^) is satisfied, then the solution 
{v, F) of the system U.l\) satisfies the following identities: 

V ■ = 0, and V ■ E^ = 0, (2.6) 

for all time t > 0. 

Proof. Following [25l |27] , we transpose the third equation of (11.11) and then apply the 
divergence operator to the resulting equation to yield 

= VjVkVjFki + VkVjVjFki. 
Using the first equation of system (II. ip . we obtain 

{VjFji)t + V ■ V(V,F,,) = 0. 
Thus, the proof Lemma 12.21 is completed. □ 

Remark 2.3. We can derive the general form of the above identity V ■ F'^ = from 
the definition of F{t,X) in (12.11) . However, the proof of the above lemma gives the 
consistency of the system. The two algebraic identities dxi ^^^^ = and = 

^ oFij Or 

{det F)F~'^ give the result of 

ax,(detFF,-^)=0. (2.7) 
Hence we obtain the following constraint on the deformation tensor F: 

^^■^d^^^] = F,r/(t,x)9xJ^^^(t,X(t,x))] 

^ det FF,-^ (t, X(t, x))dx, [-^F^j (t, ^(t, x))] 



(detF) ' ' ' 'detF 

-^dx,[Fr,^{t,X{t,x))F^^{t,X{t,x))] = 
det F 
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The key ingredient of the later proof in this paper is contained in the following 
Lemma. It shows that V x £^ is of higher order. 

Lemma 2.4. Assume that the last equality of U.3\) is satisfied and {v, F) is the solution 
of system U.l\) . Then the following identity 

^mEij — VjEim = EijViEirn — EimVlEij, (2.8) 

holds for all time t > 0. 

Proof. To prove the lemma, we will establish the evolution equation for the quantity 

^niEij — VjEim — EijVlEim + EimVlEij. 

First, by the third equation of flLlI) . we can get 

dtV^Eij + vVV^Eij + VmV-VEij (2.9) 

= Vm^kViEkj + VkViVmEkj + VmVjfj. 

Thus, we have 

dt {VmEij - VjEim) + V ■ V{VmEij - V jEiJ) (2.10) 

+ {VmV ■ VE,j - VjV ■ VE,m) 
= iy kViEkj — V jV kViEkm) 

+ ^ kViiy mEkj — VjEkm)- 

On the other hand, combining (12.91) and the third equation of (11.11) . we have 

dt{Ei^ViEij) + V ■ V{Ei^ViEij) 

= Vi [VkViEkjEim + VjViEim + VmVlEij] . 

Thus, we get 

dt{Ei^ViEij - Ei^iE.m) + V ■ V{Ei^ViEij - EijViEiJ) (2.11) 

= Vi [VkVi{EkjEim — EkmElj) + (y jViEim — V mViEij) 
+ [VmVlEij - VjViEi.m)\ . 

Combing (12101) and fl2ATD . we obtain 

dtiVmEij — V jEim + ElmVlEij — EijVlEim) 

+ V ■ V{VmEij — VjEim + Elm'VlEij — EijVlEim) 
= "^l^ kVi{EkjEim — EkmElj)\ + V kViiy mEkj — ^jEkm) 

Using the first equation of (II. ip . we have 

dtiymEij — VjEim + ElmVlEij — EijVlEim) 

+ V ■ V[VmEij — VjEim + EimVlEij — EijV [Eim) 
= Vi [VkUi{EkjEim — EkmElj) + VilVniElj — Vj-Eim)] 



On the other hand, noting (12.61) . this gives 

dtiVmEij — VjEim + EimVlEij — EijVlEim) 

+ V ■ VlVmEij — VjEim + EimVlEij — EijVlEim) 

= Vz [ - {viEkjVkEim - ViEknNkElj) + Vi{VmElj - VjEim)] 
+ V/ [Vk{viEkjEim — ViEkmElj)] 

= V; [viiVmElj — VjEim + Ek^V f^Elj — EkjV kEi^)] 
+ VlVk{viEkjEim) - VlVk{viEkmElj) 

= V/[fi(Vm-E/j — jElm + EkmSkElj — EkjV kElm)\ 

At last, by using (12. 6p and the first equation of (II. ip once again, we get the evolution 
of the concerned quantity. 

dtiymEij — VjEim + EimVlEij — EijV lEim) 

+ V ■ ViVrnEij — VjEim + Ei^VlEij — EijV lEi^) 

= W lVi(y rnElj — VjEim + EkmS kElj — EkjV kElm) 
+ ViiViEkTrNkElj — V lEkjV kElm) 

= V iViiy mElj — VjEim + Ekm^kElj — EkjV kElm) ■ 

During the calculations, we use the incompressibility conditions (12. 6p and the first 
equation of (ll.ip in the second, the third and the sixth equality. The last equality 
proves the lemma, since the above quantity will maintain zero all the time with zero 
initial condition. □ 



Remark 2.5. In order to demonstrate the mechanical background of the above lemma, 
we again go back to the definition of F{t,X) in (12. ip . Formally, the fact that the La- 
grangian derivatives commute yields the fact that dx^Eij = dx^Fik, which is equivalent 
to FikViFij{t,x{t,X)) = Fij\/iFikit,x{t,X)). Thus, one has 

FikViFi,{t,x) = FijViFik{t,x) 

which means that 

Vfc-Ejj + EikViEij = VjEik + EijWiEikit, x) 

This is exactly the result in the above lemma. However, the validity of the statement 
for any solution of the system (ll.ip is the merit of the above lemma. 

Finally, we make some simplifications for system (II. ID . 

In addition to their definitions as the elementary symmetric functions of the eigen- 
values, the invariants 7(A) of any 3x3 matrix A are conveniently expressed as 

71(A) = trA, 72(A) = ^ [{trAf - trA^\ , 73(A) = det A. 
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On the other hand, one can easily get the identity 

73(A + /) = ! + 7i(A) + 72(A) + 73(A). 

Combining the above identity with (12.51) . one can obtain the incompressible con- 
straint on E as 

trE = - det E - 72 (^) . (2.12) 
By a similar process, the incompressible constraint on E in 2-dimension case takes 

trE = -detE. (2.13) 

Next, we will consider the isotropic strain energy function W{F). We let fi{E), 
f2{E) and fsiE) represent any generic terms of degree two or higher at the origin. 
In the isotropic case, W depends on F through the principal invariants of the strain 
matrix FF"^ [llj. Define the linearized elasticity tensor as 



Suppose that the strain energy function W{F) is isotropic and frame indifferent, the 
strong Legendre-Hadamard ellipticity condition imposed upon the linearized elasticity 
tensor (I2.14p takes the form of: 

= (a^ - 2p^)5a5j^ + p\5i^5ji + 5^/^^), with a > /3 > 0, (2.15) 

where the positive parameters a and (3 depend only on W . They represent the speeds 
of propagation of pressure and shear waves, respectively. By (12.61) . (I2.14p and (12.151) . 
we have 



-dW{F) ^j,^ , d ^Wjl) 
^udFjjYi 

= («2 _ 2(3^)V,trE + /52(V ■ + V ■ E)^ + Vif2{E)a 
= P\W ■E\ + Wih{E)a, 

where we also used the assumptions that the reference configuration is a stress-free 
state: 

Without loss of generality, we assume that the constant (3=1. In particular, in 
what follows, we only consider the case of the Hookean elastic materials: V ■ fsiE) = 
V ■ {EE'^). The system is 

V ■v = 0, 

vi + v- Vv' + Vip = fiAv' + EjkVjEik + VjEij , (2.18) 
Et + v-VE = VvE + Vv. 



All the following proofs and results are also valid for general isotropic elastic energy 
functions satisfying the strong Legendre-Hadamard ellipticity condition, as those in 

3 Local Existence 

Although the proof of the following local existence theorem is lengthy, the idea is 
straight forward and had been carried out in the case of 2-D Hookean elasticity in [25] . 
For a self-contained presentation, we will carry out the similar proofs into our general 
cases. 

Theorem 3.1. Let k > 2 be a positive integer, and Vq, Eq G H''{fl) which satisfy 
the incompressible constraint U.3\) . Suppose that the isotropic elastic energy function 
satisfies the constitutive assumption Ii2.15\) . Then there exists a positive time T , which 
depends only on ||fo||_f/2 and \\Eq\\h'2, such that the initial value problem or the periodic 
initial-boundary value problem for ( f i. ij) (or ^2.1^) ) has a unique classical solution in 
the time interval [0,T) which satisfies 



for all j, a satisfying 2j + \a\ < k. Moreover, if T* < +oo is the lifespan of the 
solution, then 



Proof. By the Galerkin's method originally for standard Navier-Stokes equation [37] 
and later modified for different coupling system [23], we can construct the approximate 
solutions to the momentum equation of v, and then substitute this approximate v into 
the transport equation to get the appropriate solutions of E. To prove the convergence 
of the sequence consisting of the approximate solutions, we need only a priori estimates 
for them. For simplicity, we will establish a priori estimates for the smooth solutions of 
fl2.18p . Therefore, let us assume in the rest of this section that {v, E) is a local smooth 
solution to system fl2.18p on some time interval [0,T). 

In this paper, || • || will denote the L^(f2) norm, where C will be either an 
n-dimensional torus T", or the entire space i?" for n = 2 or 3, and (■, ■) the inner 
product of standard space L^(f2)'^ with d G {1, 2, 3, 4, 9}. 

The original system fll.ip possesses the following energy law: 



(3.1) 




(3.2) 




vf+ [ {W{F) -W{I))dx) + f^WVvf = 0. 



(3.3) 



Equivalently, for (12.181) . the corresponding energy law will be: 




(3.4) 
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which follows from the third equation of (12.181) and the incompressibility. 

The following well-known interpolation inequalities are results of the Sobolev em- 
bedding theorems and scaling techniques [2], [23] . They will be frequently used in the 
following higher order energy estimates. 



Lemma 3.2. Assume v G W^''^{Q), k > 3. The following interpolation inequalities 
hold. 



1. Forl<s< k, 



2. For2<s< k, 



II^IU* 

\Vv\\Li 
\Vv\\Li 

\Av\\l4 



\\v\\l°° 
||f ||l°° 

I Vf II LOO 

I Vf IIloo 



< C'llt; 

< C\\v 

< C\\v 

< C\\v 

< C\\v 

< C\\v 



i-^||V"t;||A, nCR^, 
i-ii||V^t;||3^, nCR^, 

i-4(ii2) ||v*At;||4{^, Q C R^ 



< c\\v\\^--^\\\/'v\\-^, ^]ci^2^ 

< C\\v\\^-^\\V'v\\^, QCR\ 

< C||t>f-^||V"Vt;||^, fiCi?2, 



< C||i;||^"^(^||V"Vt;||^(^, ncR^. 
The following two propositions can be found in [21 [161 123] . 



Proposition 3.3. If g : i?" — > R is a smooth function with g{0) = 0, then, for any 
positive constant k, g{v) G fl H'^ if v e L°° f] and 

\\g{v)\\Hk < C\\v\\Hk 

for some constant C depending only on g, k and ||t>||Loo. 

Remark 3.4. The above proposition is only used in the cases of general elastic energy 
functions. 



Proposition 3.5. Assume that f,g^ H^{Q). Then for any multi-index a, \a\ < s, 
we have 



|V"(/^7)|| < c(||/|Uoo||v^^7ll + ll^?IU==l|V7||), 
|V°(/(7) - fV^gW < C^(l|V/|U^||V^-i^7ll + I|V(7||l-||V^-7II). 

for some constant C depending only on n. 



(3.5) 
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We divide the proof of the theorem 13.11 into two parts. 

Step 1. estimate. 

Integrate (13.41) over (0,t), one obtains 

(||t;f + pf ) + 2/i f \\Vv\\^ dt = (||t;o|r + ll^oH')- (3.6) 

By taking the LP' inner product of the second equation in (12.181) with Af , using 
Lemma 13.21 and integration by parts, we have 

^^\\^vr 

= {vt, Av) + {v ■ Vv, Av) + {Vp, Av) 
- {EjkVjEik, Av') - (V ■ Av) 

< C\\Av\\ (\\vt\\ + ||^^||L-||Vt;|| + ||E||loo||V^|| + ||V^||) 
<C||Ai;||(||i;i|| + ||i;f-^(2)||^^||^?(2)||^^|| 

+ (||Ef-^(2)||AE||''(2) + 1) (ll^ll + IIA^II)) 
<^^^\\Avr + g{\\vtl\\Vvl\\AE\\). 

where 6{s) (0 < 6{s) < 1) represents a generic function which is determined by Lemma 
13. 2[ and g{-,-,-,) represents any generic nonnegative and increasing function of its 
variables. Thus, we have 

\\Avf <g{\\vt\\,\\Vv\\,\\AE\\). (3.7) 

In the meantime, by taking the inner product of the second equation in (I2.18P 
with Vf, using Lemma [3.21 and integration by parts, we obtain 



= -{v Vv, vt) - (Vp, vt) + {E.kVjEik, vl) + (y-E, Vt) 
< \\\/vt\\{\\E\\l, + \\v\\% + \\E\\) - {E.kVjEjk^vD 

<^\\Vvtr + g{\\vtl\\VvUAE\\). 

In order to get the first inequality of the above computation, we used the constraint 
on E, which is due to the incompressibility, in Lemma [2. 2[ 

Next, taking t derivative of the second equation in (12.181) . and then taking the LP 
inner product of the resulting equation with vt-, we can apply Lemma [3^ and integration 
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by parts to obtain 
1 d 



= -{dt{v ■ Vv),vt) - {Vpt,vt) + {dt{EkjdkE,,),vl) + {dtdjE,„vl) 
= -{vt ■ Vv, vt) - {dt{EkjEij), dkvi) - {dtEij, djvl) 
= iv(g) Vt, \/vt) - (dtiEE^), \/vt) - (dtE, \/vt) 

<l|Vi;t||(||t;t||||t;||L^ + ||i?*||ll^llL- + ll^t||) 
<l|V^;t||(||t;t||||^;|r-'(^)||At;r(2) 

+ \\EMEr-'^'^\\AE\r^ + \\Et\\) 
<^\\Vvtr + g{\\vt\\,\\AE\\,\\Et\\,\\Av\\). 

On the other hand, from the transport equation of fl2.18p we have 

ll^tll < ||^||L-||Vt;|| + ||t;|Uoc||VE|| + ||Vt;|| 
< C\\E\\^-'^^^\\AE\f^^^\\Vv\\ + ||Vt;|| 

<^(||V^;||,||AE||,||A.;||). 



Substitute (13. 7p into the above inequahty, one has 

\m<g{\\vt\\,\\Vv\\,\\AE\\). (3.10) 
Plugging (13.71) and (13.101) into (13. 9p . one arrives at 

~\\vtr + YllVi^tf < ^7(||'..||, ||V'.||, \\AE\\). (3.11) 

Noting (13.81) and (13.111) . it is clear that the key now is the estimate of the term 
||A£'||. It follows from the transport equation in (I2.18P that 

i^llABf (3.12) 

= -{A{v ■ VE), AE) + {A{\/vE), AE) + {VAv, AE) 

= - {A{v ■ VE) - V ■ VAE, AE) + {A{VvE), AE) + {VAv, AE) 

< C||AE||(||A^||||V^|Uoo + \\Av\\l4VE\\l. + \\VAv\\\\E\\l^ + \\VAv\ 

< C'||A^||(||AE||||VAt;f(2)||^||i-e(2)^ ||y^^||||^^||e(2)||^||i-^^^^ 

+ ||VAt;f (i)||t;||i-^(i)||A^f + ||VAi;| 

<^?(||AE||)||VA^||. 
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On the other hand, by applying V to the momentum equation in (12.181) and then 
taking the inner product of the resulting equation with VAf , we can get 

/x||VAt;f 

= {Vvt, VAv) + {V{v ■ Vv), VAv) + (VVp, VAv) 
- {W{EkjdkE,,),WAv') - {\/djE,,,VAv') 

< C\\VAv\\[\\Vvt\\ + \\Av\\\\v\\l^ + ||Vt;||i4 
+ \\AE\\\\EU^ + \\VE\\% + \\AE\\^ 

< C\\VAv\\ (^WVvtW + \\Av\\\\Av\f^^^\\v\\^-'^^^ + ||t;f At;||2^(^) 
+ ||AE||||AEf (2)||^||i-e(2) ^ ||^||2-2e(i)||^^||2e(i) ^ ||^^||^ 



Using (13. 7p again, it yields 

||VAt;|| < C\\Vv,\\ +9{\\v,\\, \\Vvl \\AE\\). (3.13) 
Insert (I3.13P into ( 13.12p . one concludes that 

Ij^W^Ef < ^\\Vv,r + g{\\v,l \\Vvl \\AE\\). (3.14) 
Combining (13. 8p . (13.111) with (13.141) . we arrive at 

|(iiAEip + f^\\vvr+\\v,r) + {f,\\vv,r+\\v,r) (3.15) 

< ^(||.;,||,||V.;||,||AE||). 

It follows from the momentum equation in (12.180 that 

\\vt{0,x)\\<C{\\vo\\m,\\Eo\\m). (3.16) 

(I3.15p . (I3.16P and the Gronwall's inequality guarantee the fact that there exist positive 
constants T,Mq, depending only on ||t>o||/f2, ||£'o||i:/2 such that 

{WAEf + fi\\Vvf + Wvtf) + r ifiWVvtf + Wvtf) ds < Mo. (3.17) 

Jo 

Returning to (13.70 and (I3.10p . we find that 

\\Av\\<giMo), m\<g{Mo). (3.18) 
And recall (13.131) . we can obtain from (I3.17P that 

rp 

I \\VAvf ds<g{Mo). (3.19) 
Jo 
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By (13.61) . fl3.17p - fl3.19p . we conclude that there exists a sufficiently large positive con- 
stant M depending only on ||i;o||_ff2, ||ii^o||//2 such that 



T 



{Ml^ + ml. + \\vtr + \\E,r) + / (||V^;||^. + ds < M. (3.20) 

We complete the proof of (13.11) when k = 2. 

To prove (13. 2p . we assume that T* < cxo is the maximal existence time and 

T* 

\\VvfH2 dt < +00. (3.21) 

Go back to (13.120 . we can use Gronwall's inequality to get 

||AE||<+oo, 0<^<T^ (3.22) 
On the other hand, by fl3.22p and the transport equation of E in (12.180 . we have 

\\Etf < {\\E\\lo.\\Vv\\ + \\v\\l^\\VE\\ + \\Vv\\Y 

<K+\\Vv\\l2, 0<t<T*. 
Thus, by fl3.2ip . we obtain 

\\Etfdt <oo. (3.23) 



If we go back to (13. 8p . using (13.61) and (13.221) . we have 

IIVT l|2 



5^l|Vt;r+||t;tr (3.24) 
<C\\VvME\\'-'^'^\\AEf^'^ 
+ ||i;f-^«||Vt;f« + ll^ll) 
< ^\\Vvtf + C\\Vvf + C. 



Similarly, by and IK22^ . IK9h will give 
1 d 



2^^n^H,' + /i||Vt;*f (3.25) 
<l|V^;J(|kt||||t;|r-''(2)l|At;r(2) 

+ m\\\Er-'^'^\AEr'^ + m\) 
<^\\vv,r+c\\v,r{\\Avr+i) + \\E,r. 

Combining fl3.24p and fl3.25p . we have 

V^f + \\v,r) + \\v,r + f^WVv.r < C||^;,f (IIA^;^ + 1) + \\E,r. 
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With (13.211) and fl3.23p . we can use Gronwall's inequality to get 

iW^vf + \\v,f) + r {WVikf + \\vtf) dt < +CX). (3.26) 
Jo 

Insert ((322]) and (13:261) into (EZD and (13101) . we get 

||Aw|| < +00, ||Et||<+oo, 0<t<T*. (3.27) 
Combining ([SJ]), (131211) -( 131221) and (131261) -( 13:271) . we get 

(11^11^2 + 11^11^2 + WvtW^ + ll^tf ) + / (||Vu||^2 + \\vt\?H^) ds < +00, 

which contradicts with the assumption that T* is the maximal existence time, which 
in turn proves the equation (13. 2p when k = 2. 

Step 2. Higher order energy estimate. 

The proof for /c > 2 is an induction on k. Assume the theorem is valid for integer 
k. In other words, we have 

\ d^^F G (0, T; i/fc-a^-l"! (f])) . ^"^-^^^ 

for all i,a satisfying 2i + |a| < k, T being determined as that in step 1. Namely, for 
all i,a satisfying 2i + \a\ < k, we have 

{\\dy''-^'vf + \\dy'"^'Ef) + / \\dy''+^'^'vf dt < +00. (3.29) 

Jo 

Here and in what follows the summations are performed over repeated indices i regard- 
less of their position, as we had assumed before. Our goal is to prove that the above 
results are valid for all j, a satisfying 2j + |a| < + 1, which are equivalent to: 

{WdlV'+^-^'vf + WdlV'-^^-^^Ef) + r WdlV'+'-^^vf dt < +00. (3.30) 

Jo 

where the summation over j is from to | if is an even number, and from to 
if k is an odd number, respectively. 

First, we assume that k is an even number and (13.291) is satisfied. By applying 
^jyfc-2j ^Yie second equation in (12.181) . we have 

dlV^-^'vt + dlV^-^^{v ■ Vv) + diV^-^^p (3.31) 
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By taking the LP' inner product of the equation (13.311) with dl'^^W^ ^•'f , < j < | 
and using integration by parts, we get 

|^||a^V'=+^-%f + ||a^V^-2%f (3.32) 

< \\di^^V''~^'v\\ (^\\dlV''-^\v ■ Vv)\\ 
+ \\d{v''-^^V ■ {EE^)\\ + ||a/v'--'^V ■ E\\y 

Applying Lemma [3^ the induction assumption 03.290 yields 

\\dlV'-'\vVv)\\ 

= \\V\v-Vv)\\ + \\dhv^v)\\+ WdiV'-^'iv-Vv^ 

0<j<l 

<\\v\\Lo^\\V''^^v\\ + \\v\\L^\\dfvv\\ + C J2 \\^^v\\l4^''^^'^v\\l^ 

0<«<f 

o</<| 

0<j<^,(l,n)j^U,k-2j) 

Further computation shows that: 

\\d{v'-'^ivVv)\\ (3.33) 

<c{i + ii9^v'=+i-2i^ii) +c Yl ii9rv^-2i-n^||i-e(i) 

0<i<|,(i,n)7^{j,fc-2j) 

<C(l + ||9/v'=+^-%||). 
In a similar way, we also have 

||^,yfc-2iy . {EE^)\\ + \\d{V^-^^V ■ E\\ < C{1 + ||^V'+^"'^'E||). (3.34) 
Putting these estimates fl3.33p - fl3.34p into fl3.32p . we obtain 

i|a/V'=+^-%f + Wdt'V'-'^vf (3.35) 

(JjL 

< C{1 + ||^V''+^-'-''t;||' + Wd'tV'+'-^'EWY. 
17 



By taking the inner product of the equation fl3.3ip with dlV^ ^•' Af , for < j < 
I and using integration by parts, we get 

+ WdiV'^'^^V ■ {EE^)\\ + \\div''-^'v -Ell). 
Noting ^Ml, dSJl, we get 

Applying 5^ V^"'"^"^-' to the third equation in (12.181) gives 
^jyfc+i-2i^^ + d{V''+^-^^{v ■ VE) = dlV^+^-^^{VvE) + dlS/^+^~^^S/v. (3.37) 

Now, we take the inner product of (13.371) with S/V'^+^^^-'-E', < j < |, and use 
integration by parts: 

\jt\ny'^''''Er (3.38) 

= {div''^^-^^{VvE),div''^^^^'E) + {div^^^~^h,div''^^-^^E) 

< \\d{v''+^-^'E\\{\\diV''+^-^'{VvE)\\ + WdiV^'+^-^^vW) 
+ \\diV^+^-^^E\\\\dlV^+^-^^{v ■VE)-v Vd{V''+^-^^E\\. 
By a similar process as in (I3.33p . we can have 

\\diV''+^-^'{VvE)\\ (3.39) 

< \\d}v{VvE)\\ + \\V'^\VvE)\\ + \\d{v'^'-'^v\\\\EU^ 
+ \\Vv\\Lo.\\div'+'''^E\\ 

0<j<|,(/,n)7^(0,0),(i,fc+l-2j) 

< C{1 + \\VAv\\){\\d{v''+^~^'v\\ + \\dlV''+^-^^E\\). 
On the other hand, we can estimate the last line of (I3.38P as follows 

\\dlv^+^-'''{v ■VE)-v Va^V'^+'-'^'^ll (3.40) 

< ||Vw||loo||V^V^-2^'E|| 

+ ^ E E ii5rv'=+^-2^-"t;iu4||9^v"+^E)iU4 

{l,n)^ {j,k-2j),{j,k+l-2j) 
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Combining (Km with ([S3HD-(IS3nD, we have 



j^WdiV'^'-'^Ef (3.41) 
< C{l+\\VAv\\){\\d{v''^^-^'v\\^+\\diV''+^~^'Ef) 

Combining this formula fl3.4ip with fl3.35p . we can conclude that 

< C{1 + \\VAv\\){\\d{v''+^-^'vf + II^V'^+^-'^'Ef ) + C. 
Noting that |Vf |loo dt < oo, we can apply Gronwall's inequality to get 

(i|a/V'=+^-2%f + \\d{V''-''-^'E\\^) (3.42) 

+ / ||^'^V=-%f ds<M, 0<t<T 
Jo 

where M depending only on ||fo||j/fc+i and H-Eolli/fc+i- Moreover, by (13.360 . we have 

r Wdiv'^'-^-^^vf ds < M. (3.43) 
Jo 

Together, (13.421) and fl3.43p imply fl3.30p when k is an even number. 

We now assume that k is an odd number and k > 3. Applying dj'V'^~^^~^^ to the 
second and third equation of fl2.18p . we have 

Now we take the inner product of the first equation in the system fl3.44p with 
^yfc+i-2jy^ where < j < -^y^, integration by parts yields the following 

+ {d{V''+^-^^V ■ {EE^),d{V''+^-^'v) 
+ (^ V ■ E, c^^V'^+'-^^v) 

< Wdiv'+^-^^vW^WdiV^^'-^'iv^v^ + \\div'+'-^'{EE'')\\ 
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where we used Lemma 13.21 and the induction assumption. In summary, we have 

V^+^-^f + fiMV'+'-'^vf (3.45) 

at 

Similarly, we will take the inner product of second equation of fl3.44p with 
^jyfc+i-2j^^ g < j < ttl and use integration by parts. The similar derivations as in 
(CT?]) will give us 

< \\df>W^+^~^' E\\ {v ■ WE) - V ■ dlW^+^-^^WE\\ 
M\diV^+^"^^{VvE)\\ + \\div^^^^^^v\\) 

< C||^V'+^-'^'E|| (l + (1 + \ VE\L^)\\d{v''^^-^^v\\ 
+ (1 + \Vv\Lo.)\\d{V''+^-^'E\\ + \\d{V''+^-^^v\ 

Employing the induction assumption, we have 
1 d 



< c{i + ||^v^+^-%f + \\d{v'+'~''E\\') + ^||a^v'=+2-%f . 

Combining (13.451) with (13.461) . we obtain: 

< c(i + Wdlv^+^'^^vf + \\dlv^+^-^' Efy 

Again, we use the Gronwall's inequality to deduce that 

+ ^ \\dlV''+^-'^^vf ds<M, 0<t<T. 

This conclude the proof of (I3.30p when k is odd. 

Putting all these results (13.421) . (I3.43P and (I3.47p together, we have proved (13. ip 
and completed the proof of Theorem 13. 1[ 

□ 
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4 Global Existence 



We now turn our attention to the proof of the global existence of classical solution for 
system f l2.18p . A weak dissipation on the deformation F is found by introducing an 
auxiliary function w below. The way of defining such a function reveals the intrinsic 
dissipative nature of the system. 

To avoid complications at the boundary, we only present the periodic case Vl = T" 
and the whole space case Q = i?". In fact, the case of smooth bounded domain can also 
be treated at a more lengthy, but no more difficult procedure than the proofs presented 
here. 

Unlike those previous results in viscoelastic literature [261 [291 [30] , the main difficulty 
lies in the apparent partial dissipation structure of the system fl2.18p . 

On the other hand, it also lacks the property of scaling invariance. The presence of 
viscosity on v gives a big obstacle to utilize the combination of Klainerman's generalized 
energy estimates and weighted estimates [ISl [13 |321 1321 [M]- 

The main contribution of our work is to reveal the fact that the incompressibility of 
system f l2.18p will provide us enough information for the proof of the near-equilibrium 
global existence of classical solutions. 

In the 3-D cases, the term V x is in fact a high order term! We recover the results 
obtained in p5|, where we avoided using this fact by the introduction of the auxiliary 
vector and then det V0 = 1 is enough to prove the near-equilibrium global existence 
of classical solutions in 2-D case. 

We start the proof by applying A to the transport equation in (12.180 and then 
taking the inner product of the resulting equation with AE, 



Next we apply A to the momentum equation in (12.181) and then take the inner 




(4.1) 
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of the resulting equation with At> to deduce that 

-^\\Avf + fi\\VAvf (4.2) 

= -{A{v ■ Vv),Av) + (AV ■ {EE^),Av) + (AV ■ E,Av) 

< C\\Av\\\\Av\\\Vv\lo^ + C\\AE\\\E\l^\\VAv\\ - {AE,VAv) 

< C{\\v\\h2 + \\E\\h2) (||Vi;f + \\VAvf + \\AE\ 
-{AE,VAv), 

where in the first inequahty, we used Proposition I3.5[ 
Combining (14. ip with (14.21) . we arrive at 



1 d 
2di 



(^11 At;f + II AEf ) + fi\\VAvf (4.3) 



< C{\\v\\h2 + \\E\\h2) (||Vi;f + ||VAt;f + ||AE| 



In order to extract the dissipative nature of the system, we want to combine the 
hnear terms on the right hand side of the momentum equation in (I2.18p . We introduce 
the auxihary variable w as follows: 

w = Av + -V ■ E. (4.4) 

The system (I2.18p will give the reformed equation: 

wt + A{v ■ Vv) + -V -{v ■ VE) + AVp (4.5) 
/i 

= fxAw + AV ■ (EE^) + -V ■ (yvE) + -Av. 

By taking the inner product of the resulting equation with w, we find 

^|lkf + /^l|V«;|r (4.6) 
= -{Aiv ■ Vv) + -V -iv- VE),w) 

- (AVp,w) + -(Av,w) 
/i 

+ -{V ■ {VvE),w)) + {AV ■ iEE^),w). 
A* 

Now let us estimate the right side of (14. 6 p term by term. First of all, the first term 
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can be estimated as 



I - (A(f ■ Vv) + -V ■ (f ■ VE),w) I 

< I (f • Vw, w) I + I (A(t> • Vf ) — V ■ VAt>, w) 



+ -\{V-{v- VE) -v-VV-E,w 

< I {A{v iS)v) -vS) Av, Vw) I 
+ ■ (v ® E) - V (^V ■ E,Vw 

< c(^\\\/v\\l, + ||At;||||t;||L.o + ^\\\/v\\\\E\\lo.^\\Vw\\ 
<C{1 + ^) HMh'^ + \\E\\h.) {llVwf + \\Vvr + \\VAt 

Next, we estimate the last term on the right hand side of (14.61) as follows: 



- (V ■ {VvE),w)) + (AV ■ {EE^),w) 
/i 

<C\\Vw\\\\Eh^{\\AE\\ + kvv 



<C{l + -)\\EyJ\\V 



'w\ 



|Vt;f + ||AE|p 



Here we used Proposition [ST 
It is rather easy to get 



(4.7) 



(4i 



(4.9) 



At last, let us estimate the term (VAp, w). Noting that V ■ v = and (12. 6p . by 
applying the divergence operator to the momentum equation of (12.181) . we get 



By Lemma [3.21 we have 



\\VE\\l < 



Ap = VjEikViEjk - ViVjVjVi 



||E||^||AE||i<||E||^.||A£;||, in 



\\E\\-qAE\\i < \\E\\h4^E\\ in 
This gives us the following estimates: 



\{AVp,w)\ < \\Vw\\[\\VE\\% + \\Vv\\l,) (4.10) 
< C{\\E\\h2 + \\v\\h2) {\\Vwf + WVAvf + \\Vvf + \\AEfy 
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Combining all the above tedious but standard estimates fl4.6p -( HrTU]) together, we 
arrive at the following important energy inequality for the auxiliary variable w. 

^||wf + /i||Vwf (4.11) 

<^jvvr + c{i + ^){\\vy. + \\Ey.) 



X (^||AEf + II Vt;f + II VAf f + ||Vt 

The key here is to estimate the tern AE. Recall the Hodge decomposition 

= VV ■ ^ - V X V X 
Here is the place that we will use (12.81) and (14. 4p to obtain the following estimate: 

\\AEf = II VV ■ Ef + II V X V X Ef (4.12) 

< 2fx^(^\\Vwf + ||VAt;f ) + ||V X V X Ef 

< 2fi^(^\\Vwf + WVAvf'^ + C\\E\\lj4AEf, 
which gives us the bound 

WAEf < Cfi''(\\Vwf + ||VAt;f ) (4.13) 

provided H-EHi^z < 

With the above result, we are ready to employ the same method as that in [25] 
to prove the global existence results. Combining (14.31) . (14.111) with (14.131) . we finally 
arrive at 



(^\\wf + \\AEf + II At;f ) + ^(^11 Vwf + WVAvf'^ 



dt 

< C{i? + ^) {\\v\\h2 + \\E\\h2) (II V^i;f + ||Vt;f + || VAt;f ) + ^11 Vt;||l 
Thus, if the initial data is sufficiently small, we can find some T* > 0, such that 

for all < t < T*. Moreover, in this case, 

(^llVwf + ||VAt;f)rfr (4.15) 

1 / \ C 

<C{fi' + -,)i\\vo\\l. + \\Eo\\l.]+-^ / 11^^ 
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holds for all < t < T*. Noting the original basic energy law (13. 6p . we have 

1 

{\\EfH2 + \\v\\l,yt) + 12 ||Vt;|||2t/r<C(^2 + —)(||t;o||^2 + 11^011^2) (4.16) 
holds for all < t < T*. iKW\ (I4A6|) imply that if 

ll^oll^. + WEoWh < 8^3(1 +^6) (1 + ^3)2' (4-17) 

then (14.141) is still true with < being replaced by < for all < t < T*, which implies 
that (I4.14P is true for all the latter time with the uniform constant C independent of 
t and fi. Moreover, from ()4.16p . we have 

/"°^ u2 
11^11^2 + 11^11^2+/^/ \\^v\\hdt< 



2C(/i3 + l)- 

This together with the local theorem 13.11 gives the following global existence of near- 
equilibrium classical solutions for system (12.180 . 

Finally, we state the theorem in the lightly more general cases. The proof is exactly 
the same as the case of (I2.18p . 

Theorem 4.1. Consider the viscoelastic model U.l\) with the initial data / li.^j) in the 

whole space or n-dimensional torus T", for ri = 2,3. Suppose that the initial data 
satisfies the incompressible constraint and the strain energy function satisfies the 

strong Legendre-Hadamard ellipticity condition \2.15\) and the reference configuration 
stress free condition ( [^. iTp . Then there exists a unique global classical solution for 
system which satisfies 



POO 

\E\\h2 + \\v\\jj2 + \\\/v\\jj2dt < 

Jo 



2C(/x3 + 1) 

if the initial data Vq, Eq G H^{Q) and satisfies the condition: 

II l|2 I II 7? I|2 ^ r'' 



M(l +/il2)' 

where k is an integer and k >2, M > 8C^ is a large enough constant. 



5 Incompressible Limits 

In numerical simulations and physical applications, one often views the incompressible 
system as an approximation of the compressible equations when the Mach number is 
small enough. Thus, it is of interest to see whether the solution to the incompressible 
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system can be obtained as the incompressible limit of the corresponding compressible 
system. Moreover, incompressible limit is also very important in the mathematical 
understanding of different hydrodynamical systems and has been extensively studied 

[MlEolEaiSlj. 

The corresponding compressible viscoelastic system takes of the following form: 
dtp + V ■ Vp + pV ■ f = 0, 

dtv + vVv + X^^^V p = ^^{Av + ■ v)) + \V ■ {pFF^), (5.1) 
dtF + vVF = VuF. 

where p{p) is a given equation of state independent of the large parameter A with 
p'{p) > for p > 0, and A the reciprocal of the Mach number M. For simplicity, we 
only concern the Cauchy problem of system (15.11) . The initial data takes 

p^(0,x) = l+p^(x), t;^(0,x) =t;o(x)+^Jo^(x), F^(0, x) = Fo(x) + Fo^(x). (5.2) 

where p'^(0,a;), F'^(0,x) satisfy 

p^(0,x) detF^(0,a;) = 1, 

Vq{x), Fq{x) satisfy the incompressible constraints (11.31) and ])q{x), Vq{x)^ ^ire 
assumed to satisfy 

\\pI{x)\\s<5,/\\ \\vl{x)\\s+i<5o/\ \\F^{x)\\s<5o/\. (5.3) 

Here 5q is a small positive constant and s is an integer with s > 4. 
For the above system, we can state the following theorem: 

Theorem 5.1. The global classical solution for system U.1\) - I[I7^) can be viewed as the 
incompressible limit of system l{5.1\ )- l[5^) if U.3\) . I{2.15\) . ^2.11^ and h5. 3|) hold and 
the incompressible initial data satisfies 

ll'^ollffs + II-E'oIIh'' — ^0 

for a sufficiently small constant Eq. 

The proof of Theorem 15.11 relies on the following Lemma 15. 2[ namely, the uniform 
energy estimates with respect to the parameter A, which was proved in [22] in 2-D case. 
The methods to prove the lemma, as well as the theorem, are very similar in the 3-D 
cases here. We will not repeat the process and want to refer to [22] for details. 



Lemma 5.2. Consider the local solutions of the compressible viscoelastic model 1(5. 
/ TOj) under the constraints / fO|) . \2.13\) . \2.11/\j and [5^) . Then the solution {p^, v^, F^) 
to system / I5. i)) -/ f3T^) satisfies the following estimates 

E,iV\t))+pJ^\\Vv'\\Ut<Ceo, . . 

Es-iidtV\t)) + p Jl^WVdtvYs-i dt<CexpCt. ^ ' ^ 
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for any t G [0,T'^] and a universal constant C independent of A if the initial data 
satisfies 

WvoWls + \\Eo\\ls < Sq. 
Here is a small enough constant and the energy Es{y^{t)) is defined as 

E.(y^t)) = ||A(p^-l)||^. + ||.;^||^. + ||E^||^.. 

Moreover Tx oo, as \ ^ +00. 
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